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1. SUMMARY 


This final report summarizes the recent results obtained by the principal investigator 
and his coworkers on the robust stability and control of systems containing parametric 
uncertainty. The starting point is a generalization of Kharitonov’s theorem obtained 
by Chapellat and Bhattacharyya in 1989 and its generalization to the multilinear case, 
the singling out of extremal stability subsets and other ramifications now constitutes 
an extensive and coherent theory of robust parametric stability that is summarized 
in the results contained here. 

2, INTRODUCTION 

The stability of a linear time invariant continuous time feedback control system is 
characterized by the root locations of its characteristic polynomial £(s); for stability 
the polynomial S(s ) must be Hurwitz i.e. have all its roots in the open left half of the 
complex plane. Since control systems operate under large uncertainties it is important 
to determine if stability is robust, that is, preserved under various perturbations. 
Despite its practical importance the subject of robust parametric stability lay dormant 
for about 100 years since the Routh-Hurwitz criterion was developed. The field was 
revived with the advent of Kharitonov’s theorem which appeared in 1978. In 1989 
Chapellat and Bhattacharyya generalized this result to make it applicable to control 
systems. This generalization has given rise to a great many useful and insightful 
results related to stability margin calculations, mixed parametric and unstructured 
uncertainty, nonlinear perturbations mixed with parametric perturbations, gain and 
phase margin optimization, development of Bode and Nyquist envelopes and classical 
design theory for robust systems etc [1,2, 3, 4, 5, 6 , 7, 8 ]. This report gives a summary of 
these results without proofs. We expect these results to have an important influence 
on future developments in this field. 

3. LINEAR INTERVAL SYSTEMS 

We begin this section by describing the theorem of Kharitonov which deals with the 
family of real polynomials 

£(s) — & l $ + & 2 S 2 + • * 41 + ( 1 ) 

with coefficient vector 6 = [£o 3 — £ n ] lying in the box 

A = [® 0 ,»o] x [xi.yi] x ••• x [x n ,y n ], (2) 

The Kharitonov polynomials associated with the above family of interval polynomials 
are defined as 

K 1 (s ) = Xo + xis + t/ 2^ 2 + j/ 35 3 + • • • 

K 2 (s) - x Q + yxs + y 2 s 2 -t- x 3 s 3 + ••• 

A 3 (s) = Vo + xys + x 2 s 2 + y 3 s 3 + ••• (3) 

K 4 (s) = yo + /3 4 s + X 2 S 2 + X 3 S 3 + • ■ • . 


3 



Theorem 1. (Kharitonov’s Theorem [9]) The family A contains only Hurwitz 
polynomials if and only if K l (s), K 2 (s), K 3 (s) and K 4 ( s ) are Hurwitz. 

This remarkable theorem unlocked the door leading to the development of a large 
number of interesting results in the area of real parametric uncertainty. However 
Kharitonov’s theorem itself is of somewhat limited applicability in control problems. 
To explain this consider the control system shown below in Figure 1. 


+ 


a 


F(.‘) 


u 


G(s) 


y 


Figure 1. Feedback System 


Let 

F(s) 

and N(s ), D(s) and Fi(s), i = 
the above system is given by 




G(s) : — 


N(s) 


( 4 ) 


F 2 (s) D(s) 

1,2 are polynomials. The characteristic polynomial of 


8{s) = Fi(s)N(s) + F 2 (s)D(js). 


( 5 ) 


In (5) the Fi(s ) may denote fixed polynomials corresponding to the controller and 
N(s), D(s) may be uncertain polynomials corresponding to the plant. Kharitonov’s 
Theorem cannot deal with this situation because of the assumption implicit in it’s 
statement that the coefficients of 8(s) perturb independently. Motivated by these 
considerations Chapellat and Bhattacharyya [10] formulated and solved the problem 
of determining the Hurwitz stability of the family 

8{s) = F 1 (s)Pi(s) + F 2 (s)P 2 (s) + F m (s)P m (s) (6) 

where the polynomials P;(s) are interval. This form of the characteristic polynomial 
(6) occurs in m-input (m-output) single output (input) systems. We shall refer to such 
families of polynomials as linear interval systems. The uncertain polynomials may of 
course correspond to the plant (perturbations) or to the controller (design parameters 
to be chosen from prescribed intervals). In the following sections we describe various 
recent results, obtained by the authors and their coworkers on the robust stability 
of the above types of control systems for both parametric as well as unstructured 
perturbations. In particular we show the importance of certain segments where the 
extremal values of various types of stability margins occur. These line segments 
capture the most important structural information for the analysis and design of 
robust control systems. 
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4 . THE EXTREMAL SEGMENTS 


We first state the result of of Chapellat and Bhattacharyya [10], [11] for the special 
case (5). F(s ) will be assumed to be fixed and G(s) to lie in an uncertainty set 
described as follows. 

Write 


N(s) : = n p s p + n p _ 1 s p f-n-js + uo 

D(s) := d q s q -+■ d q -\s q 1 + ■ - • + d\S + do 

and let the coefficients lie in prescribed intervals 

Ti; e [n t “,nt], i € {0,1, •••,?} :=£ 

die[dr,df], i € { 0 , 1 , •••,?}:= q. 

Introduce the interval polynomial sets 

W(s) := {N{s) = n p s p + n p ^s p ~ x + ■ + n 0 : n; € [n t “,nf], i <E p} 

'D(s) := {jD(s) = d q s q + d q ~is q 1 + — -f d 0 : di £ [d± ,df], i £ g} 

and the corresponding set of interval transfer functions (or interval systems) 

G ( S ) = {^7 f ] : (N(s),D( a ))eM(. S )xV(s)}. (7) 

D(s) 

The four Kharitonov vertex polynomials associated with M(s) are 

K^(s) := tiq + n^s + n^s 2 + n£ s 3 + s 4 + s 5 + • • ■ 

K 2 (s) := 71q + s 2 + nj s 3 + s 1 + s 5 + • • • 

K 3 (s) := tIq + n~s + 71^ s 2 + n£s 3 + s 4 -f 7i,T s 5 + - • • 

K*(s) := no + ti^s + 7Z2 s 2 + 713 s 3 + n^s 4 + 7i^s 5 + • • • . 


and we write 

M') := *,!(»)}■ («) 

Similarly the four Kharitonov polynomials associated with 'D(s) are denoted K^(s), 
i — 1, 2, 3, 4 and 

M') :={*■<(*). ^3W, *■?('). *!(«)} (9) 

The four Kharitonov polynomial segments associated with Af(s) are defined as follows: 

Sj/(s) ■= 

[\K' n (s) + (1 - X)KUs) : A €[0,1], (*,i) € {(1,2), (1,3), (2, 4), (3, 4)}] (10) 

and the four Kharitonov polynomial segments associated with 'D(s) are denoted 

<Sd(s) := 

[ptfiM + (1 - li)Ki(s) : p e[0.1], (i, ;•)£{(!, 2), (1,3), (2, 4), (3, 4)}] (11) 
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Following Chapellat and Bhattacharyya [10], [11] we introduce the extremal subset 
as follows: 

{tf(s)xV(s)) E = {{N(s),D(s)) : N(s) e D(s) e S v (s) 

or N(s) e Sjs(s), D(s) € M*)}- (12) 

The extremal subset of the family of interval systems G(s) is naturally defined 
as: 

g eM:={^ : (JVM, DM) € (M(s)xV(s)) e }. (13) 

D(s) 

Theorem 2. The control system of Figure 1 is stable for all G(s) £ G(s) if and 
only if it is stable for all G(s) £ G E (*s). 

We note that each element of G E (s) is a one parameter family of transfer func- 
tions and there are at most 32 such distinct elements. The above Theorem therefore 
gives a constructive solution to the problem of checking robust stability by reducing it 
to a set of (at most) 32 root locus problems. An alternative way to check the stability 
of linear interval systems would be to use the Edge Theorem [12] which requires that 
the exposed edges of the polytope of uncertain polynomials in the coefficient space 
be checked for stability. While this result also leads to a set of line segments to be 
checked the number of segments (exposed edges) increases exponentially with the di- 
mension of the uncertain parameter. In the theorem [10] the number of segments to 
be checked is independent of the dimension of the uncertainty set Af(s)xV(s). More 
importantly perhaps these extremal segments enjoy many extremal properties that 
are critical in determining stability margins. We now give an example to illustrate 
the Theorem [10]. 

Example 1. Let us consider the following single input single output plant 

n(s) _ s 3 + as 2 — 2s+/3 
d(s) s 4 + 2s 3 — s 2 + 7 $ + 1 ’ 

with 

a £ [-1,-2], (3e [0.5,1], 7 £ [0? 1]« 

There are two Kharitonov polynomials associated with n(s), namely 

0.5 — 2s — s 2 + <s 3 and 1 — 2s — 2s 2 -f s 3 
and also two Kharitonov polynomials associated with d(s) 

1 — s 2 + 2s 3 4- s 4 and 1 + s — s 2 + 2s 3 + s 4 . 

In order to check that a given controller F(s) stabilizes the entire family of plants, 
we only need to check that the controller stabilizes the following four plant segments 
making up the set G E (s): 

0.5(1 + A) — 2s — (1 + 2A)s 2 + s 3 
1 - s 2 + 2s 3 + s 4 
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A <=[0,1] 


0.5(1 + A) - 2s + (1 + 2A)s 2 + a 3 
1 + s — s 2 + 2s 3 + s 4 
0.5 - 25 


S 2 +S 3 


l + As — s 2 + 2s 3 + s 4 
1 — 2s — 2s 2 + s 3 


A € [0, 1] 


A €[0,1]. 


1 + As — s 2 + 2s 3 + s 4 
On the other hand if one uses the Edge theorem, it is necessary to check the 12 plant 
segments corresponding to 


a = - 2 , /3 = 0.5, 7 € [0, 1 ] 

a = - 2, 0 = 1, 7 € [0, 1] 

a = -l, 0 = 0.5, 76 [0,1] 

a = -l, 0 = 1 , 7 € [ 0 , 1 ] 

a = -2, 0 € [0.5,1], 7 = 0 
<* = - 2 , 0 € [0.5,1], 7 = 1 
a = -l, 0 € [0.5,1], 7 = 0 
a = - 1, 0€ [0.5,1], 7 = 1 

a G [-2,-1], 0 = 0.5, 7 = 0 

a G [-2,-1], 0 = 0.5, 7 = 1 

a £ [-2,-1], 0 = 1, 7 = 0 

a € [-2,-1], 0 = 1, 7 = 1- 

The problem of checking the stability of a line segment of polynomials was solved 
by Chapellat and Bhattacharyya in [13]; this result called the Segment Lemma is 
described next. The lemma basically checks for the occurrence of a jui root along a 
line segment of polynomials. 

4.1. Stability of Segments 

Let 5j(s), and £ 2 (s) be two polynomials of degree n and let 

Sx(s) = A*i(s) + (1-A)£ 2 (s). 

Denote the segment 

{*>(•): * €[ 0 , 1 ]} :=[«,(»),«.(.)]■ ( 14 ) 

Also write 


5i(j<*>) = i|(w)+iw^(«) 
&i{jw) = <^(w) +juf%(u;). 


where and 8°{w) are real. 
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Lemma 1. (Segment Lemma [13]) Suppose that S x (s) andS 2 (s ) are Hurwitz poly- 
nomials of degreen with positive coefficients. Then there exists an unstable polynomial 
on the line segment [5i(s), 62 ( 3 )] *ff there exists ujq > 0, such that 

6l{u>o)%(u 0 ) - - 0 


and 


< 0 

The Segment Lemma completely solves the problem of checking the Hurwitz 
stability of a line segment of polynomials. The idea of checking for a root on the 
stability boundary can of course be extended to other stability regions besides the 
left half plane. 

In general it is known that the stability of the endpoints (vertices) of a line 
segment does not guarantee that of the entire segment. It is therefore useful to know 
if there exist some simple additional conditions (simpler than the Segment Lemma) 
under which the Hurwitz stability of a segment could be guaranteed. It turns out 
that by restricting the form of the difference polynomial 

8o{») = S 2 (s) - 8 1 (s) (15) 

it is possible to conclude segment stability from vertex stability. There exist several 
known results on this problem. Peterson [14] derived the case when the difference 
polynomial is antiHurwitz (all roots in the closed right half plane), Chapellat and 
Bhattacharyya [10] dealt with the case when it is even or odd, Hollot and Yang [15! 
and Mansour and Kraus [16] proved the result for the difference polynomial being of 
the form s t (as + b)P(s) where P(s) is even or odd. The general result given below 
encompasses all the above cases; the proof may be found in Bhattacharyya [17]. 

Lemma 2. (Vertex Lemma) [17] 

a) Let £i(.s) and <$ 2 («s) he polynomials with positive coefficients and let 

8 0 (s) = A(j)5*(a5 + 6)P(s) 

where A(s) is antiHurwitz t > 0 is an integer, a, b are arbitrary real numbers, 
and P{s) is even or odd. Then stability of the segment [Si(s), £ 2 ( 3 )/ i s implied 
by that of the endpoints Si(s), S 2 (s). 

b) When Sq(s) is not of the form specified in a), stability of the endpoints is not 
sufficient to guarantee that of the segment. 

Using the above Lemma in conjunction with the linear case theorem [10] it 
is possible to show that if F{(s ) are of the same form as 6 q( 5 ) the stability of the 
extremal segments (and therefore the robust stability of a linear interval system) can 
be ascertained from the Kharitonov vertex polynomials of the system. To state this 
result let 

G*(») := { £$ ■ (,NW,D(s)) € (/C M (s)xK v (s)) } . (16) 
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Corollary 1. (Theorem 2) [17] , [1 1] Under the conditions of Theorem 2, if 

Fi(s) = A i (s)s ti (a i s + bi)Qi(s),i = 1,2 

where ^( 5 ) is antiHurwitz , ti is an arbitrary nonnegative integer, a{,b{ are arbitrary 
real numbers and is either an odd or an even polynomial, then the system of 

Figure 1 is stable for all G(s ) £ G(s) if and only if it is stable for all G(s ) £ Gk(<s). 
Moreover if the F^s) do not satisfy the above conditions stability of the system can 
not be guaranteed by verifying stability for G(s) £ Gk(s). 


This result is obviously useful in control problems where F(s) the compensator 
consists of integrators, first order lags and leads and unstable and nonminimum phase 
elements. For instance in Example 1 if we consider the problem of stabilizing the 
interval family given, with a controller of the form 


F(s) 


as “I - b 
+ d) 


it is only necessary to check that the controller simultaneously stabilizes the set of 
plants Gk(<s) shown below: 


0.5 — 2s — s 2 + s 3 

0.5 — 2s + s 2 + s 3 

1 — s 2 + 2 s 3 + s 4 

1 + s — s 2 + 2s 3 + s 4 

0.5 - 2s - s 2 + s 3 

1 - 2s - 2s 2 + s 3 

1 — s 2 + 2s 3 + s 4 

1 — s 2 + 2s 3 + s 4 

1 - 2s - 3s 2 + s 3 

1 — 2s -j- 3s 2 d" s 3 

1 — s 2 + 2s 3 -|- s 4 

1 -f s — s 2 + 2s 3 ) s 4 

0.5 - 2s - s 2 + s 3 

1 - 2s - 2s 2 + s 3 

1 + s — s 2 + 2s 3 + s 4 

1 + s — s 2 + 2s 3 -f s 4 


To state the theorem for the general case (6) we introduce some notation. For any 
positive integer n let n denote the set of integers {1,2, ■••ti}. Referring to (6) let 
d°(Pi) denote the degree of P{(s) and let p\ denote the coefficient of s l in Pi(s): 


Write p. 

—1 


P i( S ) =Pi +Pi S + 1" PW d +- pf (P,) 5 d ° (P<) - 

(p'iiPi) ' ' ' pf^ P '^) f° r eac h i 6 Hh and let 
P := \Pvl 2 


(17) 


(18) 


denote the vector of coefficients of the polynomials Pi(s ), i £ m. Each such coefficient 
belongs to a given interval: 


Pi 


e 1 = 0, e m. 


(19) 


We let Vi(s) denote the set of interval polynomials to which .P»(s) belongs and intro- 
duce the uncertainty set 


n := {(‘Pi(s)x'P 2 (s)x • • • , x7? m (s))} 
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Alternatively the uncertainty set can be described in the space of the polynomial 
coefficients. With mild abuse of notation we use II to also denote the box of uncertain 
parameters: 

n := {p I V \ G K /?'],/ = 0,.. - € m}. (20) 

Each point p G II corresponds to a particular choice of the ordered set of polynomials 
Pi(s), i G m. We write 6(.s, p) for the polynomial family (6) to display the explicit de- 
pendence of 6(.s) on p. For a given fixed set of polynomials [Fi(.s), F 2 (s), • • • , -F m (.s)] := 
F_ let A denote the family of polynomials generated by the map F_ : n => 5 as in 
(6) and obtained by letting the parameter vector p, (equivalently, the polynomials 
Pi(s)), range over the box II described in (20). In other words 

A :={*(*, p)|p€ll}. (21) 

Following the previous notation the four Kharitonov polynomials and polynomial 
segments associated with Pi(s) are denoted /C;(s) and «S;(s) respectively and these 
definitions hold for each i G m We now introduce some special subsets of II. The 
linear manifolds II/,/ G 221 are defined: 



n, := 

{(/C a (s)x • • • x£/_i(s)x<S((s)x/C I+ i(s)x- • • ,x/C m (s))} 

(22) 

Also 

A, := {6(s,p)|p G n,}. 

(23) 

Finally, let 

771 



n E := U n * 

/-I 

(24) 

and 

m 

Ae := (jAi = {^(■ s jP)Ip £ He}- 

i=i 

(25) 


Since there is a one to one correspondence between the elements of IIe and of Ae 
we refer to both sets as extremal segments. We also define the Kharitonov vertices 
K(n), of II to be the subset of all vertices of II corresponding to the Kharitonov 
polynomials of the P;(s). It is not difficult to see that the number of distinct segments 
in IIe in the most general case when all the Kharitonov polynomials associated with 
each polynomial P;(.s) are distinct, is rn4 m . With these preliminaries we are ready to 
state the main result of this section. 

Theorem 3. [10,11] F_ stabilizes II if and only if F_ stabilizes IIe- 

In the next section we show that these extremal segments are useful in deter- 
mining how close one is to instability in the parameter space over a stable set of 
parameters II. 
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5. EXTREMAL PARAMETRIC STABILITY PROPERTIES: LINEAR 
CASE 

We now turn to the question of relative stability of the family II. In other words 
given a family of polynomials II which is stable, we wish to know the “distance” to 
the closest unstable polynomial as the point p (equivalently the set of polynomials 
P f (s)) varies over the box II. Before discussing the case of a control system we deal 
with the special case of a single interval polynomial. In this case we first establish an 
important extremal property of the Kharitonov polynomials, namely that the closest 
point to instability over a stable box in coefficient space lies at one of the Kharitonov 
vertices. The proof of this result was first given in [18]. 

5.1. Extremal Property of the Kharitonov Polynomials 
Suppose that we have proved the stability of the family of polynomials 

£(s) = T Sis T 82 s 2 + • • • + 8 n s n , (26) 

with coefficients in the box 

A = [s 0 , 3/o] X [*i,»i] x x [x n ,y n \. (27) 

Write 8 = [8 0 ,8 1 ,---8 n ], and regard 8 as a point in IZ n+1 . Let ||£|| p denote the p 
norm in f? n+1 and let this be associated with 5(s). The set of polynomials which are 
unstable of degree n or of degree less than n is denoted by IA. Then the radius of the 
stability ball centered at 8 is 

p(S) = inf ||£ - u||„. (28) 

If the polynomial family is stable it is possible to associate with each element the 
largest stability hypersphere around it. We thus define a mapping from A to the set 
of all positive real numbers: 

A 7£ + \{0} 

*{•) — M 

and ask the question: Is there a point in A which is the nearest to instability? Or 
stated in terms of functions: Has the function p a minimum and is there a precise 
point in A where it is reached? The answer to that question is given in the following 
theorem proved in [18]. In the discussion to follow we drop the subscript p from the 
norm since the result holds for any norm chosen. 

Theorem 4. (Extremality property of the Kharitonov polynomials [18]). 

The function 

A 7l + \{0} 

6(s) — . p(S) 

has a minimum which is reached at one of the four Kharitonov polynomials associated 
with A. 
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The above optimal property of Kharitonov polynomials is extremely important 
and has many uses. Using it it is possible to prove the extremal property of the 
stability segments occurring in the linear case result given in the previous section. 

5.2. Extremal Parametric Property of the extremal Segments 

Consider the family A and the segments I1e and Ae which occur in the linear case 
result of the previous section. As before consider the family of polynomials lying in 
the uncertainty set II and let the coefficients of the polynomials Pi(s), p 6 i?" vary 
in the prescribed box II. Let || • || denote any norm in R n and let V u denote the set 
of points u in R n for which £(.s,u) is unstable or loses degree (relative to its generic 
degree over II. Let 

P(P) = inf ||p - u|| p (29) 

lit r ti 

denote the radius of the stability bail (measured in the norm || ■ ||) and centered at 
the point p. This number serves as the stability margin associated with the point p. 
If the box II is stable we can associate a stability margin with each point in II. A 
natural conterpart of the question posed in the previous section is: Is there a point 
in II which is closest to instability in the norm || • || and where is it? The answer to 
that question is provided in the following theorem first proved in [19]. 

As before we define a mapping from II to the set of all positive real numbers: 

n H + \{0} (30) 

p — ♦ p(p)- 

Theorem 5. (Extremal property of the Extremal Segments [19]) The func- 
tion 

n ft + \{o} 

p — ♦ p(p) 

has a minimum which is reached at a point on the extremal manifolds IIe- 

We remark that this optimality property enjoyed by the extremal segments is 
very useful to find the worst case parametric robustness margin associated with a 
given controller. 

6. PARAMETRIC AND UNSTRUCTURED PERTURBATIONS: LIN- 
EAR CASE 

We now turn our attention to problems where parametric as well as unstructured 
uncertainty is simultaneously in operation. The main results in this section will deal 
with the calculation of the stability margin for systems containing parameter 
uncertainty as defined above. In the following we will use the standard notation: 
C+ := {s £ C : Re(s ) > 0}, and i/ 00 (C + ) will represent the space of functions f(s) 
that are bounded and analytic in C+ with the standard norm, 

H/lloo = sup|/(jw)|. (31) 

ll>£R 


12 



Consider first the feedback system below where the fixed stable system with 
transfer function G(s) is perturbed by norm bounded feedback perturbations 
A G. According to the small gain theorem the perturbed system remains stable as 
long as HAGIloo < We note that can be regarded as a “complex gain 

margin” for the system. The obvious first step is to generalize this by letting G(s) lie 
in an uncertainty set. We consider the case when G(s) = belongs to an interval 
family G(,s) as in Section 4. We adopt the notation of Section III and let G/r(s) be 
the set of 16 Kharitonov systems associated with G ( 5 ), Our robust version of the 
small gain theorem can be stated as follows. 


Theorem 6. (Robust Small Gain Theorem [20]) Given the interval family G(s) 
of stable proper systems, the closed-loop system in Figure 2 remains stable for all stable 
perturbation A P such that HAPU^ < a if and only if, 


a < 


1 

max G(i)eGK (.) ||G|| OO 


(32) 


The proof of this theorem given in [20] was based on the following fundamental 
characterization of H ^ norms in terms of Hurwitz stability of polynomials. This 
lemma was also proved in [20]. 


Lemma 3. [20] Let h(s) = n(s)/d(s) be a proper (real or complex) rational function 
in Hoo(C+), with deg(d(s)) — q. Then ||fe||oo < 1 if and only if 

al) |n,| <\d q \, 

bl) d(s) + e^ 8 n(s) is Hurwitz for all 8 in [0, 2i r). 


We now give an example to illustrate the above theorem 


Example 2. Consider the following stable family G(s) of interval systems whose 
generic element is given by, 

n 0 + ni«s + n 2 s 2 + n 3 s 3 
do ~b d\s dzs 2 d$s^ 
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with 


and 


n 0 G [1,2], 71 X G [-3,1], n 2 G [2,4], n 3 G [1,3] 


do £ [1>3], di G [2,4], d 2 G [6,7], d 3 G [1,2]. 
Gk(«s) consists of the following 16 rational functions, 


9i(e) = 

03^) = 

55(5) = 
97{s) = 
9d(e) = 
9u(e) = 
9is(e) = 
9is(e) = 


1 — 3^ + 4s 2 + 35 s 

1 + 25 +Y7 2 + 25 3 ’ 

1 - 35 + 45 2 + 35 s 
3 + 25 + 65 2 + 25 3 ’ 
1 + 5 + 45 2 + 5 3 

1 + 25 + 75 2 + 25 3 ’ 

1 + 5 + 45 2 + 5 3 

3 + 25 + 65 2 + 25 3 ’ 

2 - 35 + 25 2 + 35 3 

1 + 25 + 75 2 + 25 3 ’ 

2 — 35 + 25 2 + 3s 3 

3 + 25 + 65 2 + 25 3 ’ 

2 + 5 + 25 2 + 5 3 

1 + 25 + 75 2 + 25 3 ’ 
2 + 5 + 25 2 + 5 3 


1 — 35 + 45 2 + 35 3 
1 + 45 + 7 5 2 + 5 3 ’ 
1 - 35 + 45 2 + 35 3 
3 + 45 + 65 2 + s 3 ’ 

1 + 5 + 45 2 + 5 3 
1 + 45 + 7 5 2 + 5 3 ’ 

1 + 5 + 45 2 + 5 3 

3 + 45 + 65 2 + 5 3 ’ 

, , 2 - 35 + 25 2 + 35 3 

9lo{s) - J + 4s + 7s 2 + 5 3 ’ 

2 - 35 + 25 2 + 35 s 

3 + 45 + 6 5 2 + 5 3 1 
2 + 5 + 25 2 + 5 3 


92(e) = 
94(e) = 

9e(s) = 
9s(e) = 


912 (e) 
914(e) = 
5ie(^) = 


3 + 25 + 65 2 + 25 3 ’ 

The H oo norms of the above functions are given by, 


1 + 45 + 75 2 + 5 3 ’ 

2 + 5 + 25 2 + 5 3 

3 + 4s + 65 2 + 5 3 


Halloo = 2.112, 
Halloo = 1.074, 
Halloo = 3.356, 
Halloo = 2.848, 


Halloo = 3-0, || 53 ||oo = 5-002, Woo = 3.0, 

IMoo = 10, Hsrlloo = 1-710, Halloo = 1.0, 
ll^iolloo = 3.0, llffnlloo — 4.908, j|^i2 |joo “ 3.0, 
11*41100 = 2-0, ||*b||oo = 1.509, Moo = 1.0. 


Therefore, by the above theorem the entire family of systems remains stable under any 
unstructured feedback perturbations of norm less than a = ^ = 0.19992 which 

is the smallest “complex stability margin” over the given box of parameters. 


6.1. Computation of the Structured Margin 

The converse problem is: given a prescribed bound on the level of unstructured per- 
turbations that are to be tolerated determine the amount of parameter perturbations 
permissible. 

In this case one starts with a nominal stable system 


9°(e) 


n° Q + n?5 + . . . + n°s p 
d o 0 + d° lS + ... + d°s“ 


(33) 


which satisfies ||<7°||oo = ct. A bound ^ ^ is then set on the desired level of un- 

structured perturbations. It is then possible to fix the structure of the parametric 
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perturbations and to maximize a weighted l ball around the parameters of < 7 °(s). 
More precisely, one can allow the parameters n^dj of the plant to vary in intervals 
of the form 

Tii E [n° — dj E [ d j - e/Xj, d\ j + e/Xj], 

where the weights v^\Lj are fixed and non negative. For each e we get a family of 
interval systems G(e) and its associated Kharitonov systems G/f(e). The structured 
stability margin is then given by the largest e, say e moxj for which every system g(s) 
in the corresponding interval family G(e max ) satisfies ||sf||oo £ /3. 

An upper bound e* for Cmax is easily found by letting €i be the smallest number 
such that the interval family, 


(d(.s) = do + . . . + d q s q : dj E [d° — ei/Xj, d°j + e 1 /Xj]}, 

contains an unstable polynomial. This upper bound is easily calculated, using for 
example the method proposed in [20] which is straightforward. This can be used to 
initiate a bisection algorithm. The reader may consult [20] for the details and an 
example. 

We now consider the case where the fixed transfer function F(s) in Figure 1 is not 
necessarily unity. Let G(s) be a family of strictly proper interval transfer functions, 
Assume also that we have found a stabilizing controller F(s ) for the entire family. We 
therefore have a family of stable closed-loop systems and we consider unstructured 
additive perturbations as shown in Figure 3. Here also we want to determine the 
amount of unstructured perturbations that can be tolerated by this family of interval 
plants. In order to do so we have to find the maximum of the H ^ norm of the closed- 
loop transfer function F(s)(l + G(5)F(^)) 1 over all elements G(s ) in G(s). This 
result is also reported from [20]. 



Figure 3. 


Theorem 7* [20] Given an interval family G(s) of strictly proper plants and a 

stabilizing controller F(s) for G(s), the closed loop system in Figure S remains stable 
for all stable perturbations AG such that HACzHoo < <z if and only if 

1 

“ “ max G6G , ||f(.)(l + GMFMr'lU' 
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We illustrate this result with an example. 

Example 3. Consider the following family of interval plants, 

n p (s) (3s 


9pM = 


d p (s ) 1 — s + 7 s 2 + s : 


, >3 e [1,2], 7 e [3.4,5], 


Using the theorem [10] one can easily check that the controller F(s ) = stabilizes 

the entire family. The transfer function of interest is given by, 

F(s)(l + G 0 „(s)F(s)) 1 = - - 3/?5 _ l)a 3 + + \) a 8 + S 4 • 

Following Theorem 7, we have to find the maximum H <*, norm of four one-parameter 
families of rational functions, namely, 


t - a ( s ) = 
r M (s) = 
= 

r t( s ) = 


3 (1 ~ 5 + A ^ 2 + * 3 ) A € [3 4 51 

1 -)- 3s + (A — l)s 2 + (A + l)s 3 + s 4 ’ 1 ' ’ J ’ 

3(1 - s + /xs 2 + s 3 ) e r 3 4 5] 

1 + 6s + (/x - l)s 2 + (fi + l)s 3 + s 4 ’ ^ 1 ' ’ J ’ 


3(1 — s -f 3.4s 2 + s 3 ) 


1 -f 3 us + 2 As 2 + 4.4s 3 + s 
3(1 — s + 5s 2 + s 3 ) 


7 , «' 6 [ 1 , 2 ], 


, ( 6 [ 1 , 2 ]. 


1 + 3 (s + 4s 2 + 6s 3 + s 4 
Consider for example the case of r^s). We have, 

9((1 — Au; 2 ) 2 + u; 2 (l + u; 2 ) 2 ) 


Mju »)| 2 = 


(1 - (A - l)u; 2 + u; 4 ) 2 + u; 2 (3 - (A + 1 )cj 2 ) 


2 * 


Letting t = u? 2 we have to find, 

f( . 9((1 - At) 2 + t(l + Q 2 ) 

t >o S ^3,,s] J ’ “ (1 - (A - l)t + t 2 ) 2 + t(3 - (1 + A)t) 2 ' 

Differentiating with respect to A we get a supremum at, 


or 


MO = 


A 3 (t) 


-2t + 3 + + 12 t 2 + 1 

2 1 


-2t + 3 - V4t 3 + 12t 2 + 1 
2 1 


It is then easy to see that Ai(t) £ [3.4,5] iff t 6 [^1,^2] U [^3, ^4] , where, 

h ~ 0.39796, t 2 ~ 0.64139, t 3 ~ 15.51766, t A ~ 32.44715 , 
whereas, A 2 (0 £ [3.4,5] iff t 6 [is,^] where, 

t 5 ~ 0.15488, t 6 ~ 0.20095 . 
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As a result, the maximum H x norm for r^(s) is given by, 


max (|| r 3. 4 ||ooj || r 5||oo, / sup sup f(t,X 3 

\ y te[ti,t 2 ]u[t3,t4] y te[t 5 ,f«] 

where one can at once verify that, 

9(2* - 1 - + 12 i 2 + 1) 

/(*»M0) - 2*2 + 7* _ i _ (* + 1)V4< 3 + 12* 2 + 1 ’ 

and, 

9(2* - 1 + ^4*3 + 12t 2 + 1) 

/(*, A 2 (*)) - 2*2 + 7* _ ! + (* + i)^/4* 3 + 12* 2 + 1 ’ 

This maximum is then easily found to be equal to, 

max(34. 14944, 7.55235, 27.68284, 1.7028) = 34.14944 . 

Proceeding in the same way we finally get the following result, 

0e [ i “? e[3 4s ) H ^ K 1 + ^( s ) c ( 5 )) _1 |!° 0 = 34 - 14944 > 

where the maximum is in fact achieved for /? = 1, 7 = 3.4. 

The results given above lead to the following theorem. 

Theorem 8. Let G(s) be a family of interval plants of fixed degree and let a > 0 be 
given . There exists a linear time invariant controller F(s) that stabilizes G(s) and 
that satisfies , 

sup ||F(s)(l + G(s)F(s))~ 1 ||oo < a, 

Ge G 

if and only if such a controller exists for Ge- 

We believe that the above theorem sets the stage for further investigation into 
the synthesis problem by precisely specifying the role of the controller in the robust 
stability of a family of interval systems. 

In the next section we describe the connection between parametric perturbations 
and nonlinear feedback perturbations. 

7. PARAMETRIC AND NONLINEAR PERTURBATIONS: LINEAR 
CASE 

An important stability robustness problem that involves unstructured perturbations 
is the classical Lur’e problem of nonlinear control theory. This problem considers a 
fixed linear time invariant system subjected to perturbations in the form of nonlinear 
feedback gains contained in a prescribed sector. In [21] a robust version of the Lur’e 
problem was treated where structured and unstructured perturbations are simulta- 
neously present. In this formulation the fixed linear system is replaced by the more 
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realistic model of a parametrized family of plants. The “nonlinear stability margin 55 of 
the system can be determined by finding the infimum, over the parametrized family, 
of such stability sectors. From standard results on the Lur’e problem, the size of such 
a sector can be determined by finding the infimum of the real part of G(ju ;) as G(s) 
ranges over a parametrized family G(s). In [21] it was shown how the strict positive 
realness (SPR) property for a stable family of interval systems can be determined 
from the set G K (s) of the sixteen Kharitonov systems . In addition, in the presence of 
a fixed controller that stabilizes an entire family of interval systems, the SPR property 
for the family of transfer functions F(s)(l + G(s)F(s))~ 1 is determined from a the 
extremal subset of systems. These results are described in this section and the reader 
is referred to [21] for proofs. 

We begin by giving a stability characterization of the SPR property proved in 
[21]. Let G(,s) = be a real proper transfer function with no poles in the closed 
right-half plane. 

Theorem 9 . G(«s) is SPR if and only if the following three conditions are satisfied: 

a) Re G( 0) > 0, 

b) n(s) is Hurwitz stable , 

c) d{s) + jan(s) is Hurwitz stable for all a in R. 

Based on this result it is possible to formulate a robust SPR result as follows. Let 
G(s) now belong to an interval family G(s) as in Section 3 . Given a real number 7 
we ask: Under what conditions is G(s) + 7 SPR for all G(s) in G($)? The answer is 
given in the following result from [21]. 

Lemma 4 . [ 21 ] G(,s) + 7 is SPR for every element in G(s) if and only if it is SPR 
for the 16 Kharitonov systems in G^(^). 

This result leads to: 

Theorem 10 . [ 21 ] Given a proper stable family G(s) of interval plants , the 

minimum of Re(G(ju?)) over all w and over all G(s) in G(s) is achieved at one of 
the 16 Kharitonov systems in G# (s) 

We illustrate this result with an example. 

Example 4 . Consider the following stable family G(s) of interval systems whose 
generic element is given by 

_ 1 + as + /3s 2 + -s 3 
7 + Ss + es 2 + s 3 

where 

« €[1,2], >3 e [3,4], 7 €[1,2], *€[5,6], c€ [3,4]. 


18 



G k (s) consists of the following 16 rational functions 

ri(s) 

r 3 (s) 
r 5 (s) 
r 7 (s) 
r 9 (s) 
r u (s) 

r i 3 (<s) 

The corresponding minima of their respective real paxts along the imaginary axis are 
given by, 

inf Re rAjw) = 0.1385416, ini^Re r 2 {ju) = 0.1134093, 

uj^R 

inf Re r z (ju>) = 0.0764526, in{ ueR Re r 4 (jo;) = 0.0621581, 

U/Gfl 

inf Re r s (ju) = 0.1540306, ini^R Re = 0.1262789, 

inf Re r 7 (ju > ) = 0.0602399, m{ ueR Re r 8 (ja>) = 0.0563546, 

u>ER 

inf Re r 9 {ju>) = 0.3467740, inf„ 6 * Re r 10 (ju) = 0.2862616, 

uj£R 

inf Re r u (ju>) = 0.3011472, inf w6 * Re r l2 (ju) = 0.2495148, 
inf Re r 13 (ju;) = 0.3655230, inf w6fi Re r 14 (ju>) = 0.3010231, 

Cl>£R 

inf Re r\s(ju>) = 0.2706398, inf„ eR Re r 16 (ju) = 0.2345989. 

u£R 

Therefore, the entire family is SPR and the minimum is achieved at r 8 (s). 

Turning now to the Lur’e problem let us refer to Figure 4. The class of allowable 
nonlinearities is described by sector bounded functions. Specifically, the nonlinearity 
a) satisfies 

<£(£, 0) = 0 for all t > 0 
0 < cr<j){t , a) < ka 2 . 

This implies that <£(i, a) is bounded by the lines <j> = 0 and <j> = k<r . Such nonlinearities 
are said to belong to a sector [0, k\. Referring to Figure 4, we state the following well- 
known classical result on absolute stability . 


1 + 5 + 3s 2 + s 3 , __ 1 + s + 3s 2 + s 3 

= 1 + 5s + 4s 2 + 5 3 ’ r2 ^ S * ~ l+6s + 4s 2 + s 3 ' 

1 + S + 3 s 2 + S 3 , v 1 + S + 3s 2 + 5 3 

_ 2 + 5 s + 3s 2 + s 3 ’ ~ 2 + 6s + 3s 2 + s 3 ’ 

1 + s + 4s 2 + s 3 _ 1 + s + 4a 2 + -s 3 

= 1 -f 5s + 4s 2 + s 3 ’ ~ l + 6s + 4s 2 + a 3 ’ 

1 + s + 4s 2 -f s 3 _ 1 + S + is 2 + s 3 

= 2 + 5s + 3s 2 + s 3 ’ " 2 + 6s + 3s 2 + s 3 ’ 

1 + 2s + 3s 2 + s 3 , v 1 + 2s + 3s 2 + s 3 

1 + 5s + 4s 2 + S 3 ’ 10 1 + 6s + 4s 2 + S 3 1 

1 + 2 S + 3s 2 + S 3 . 1 + 25 + 35 2 + 5 3 

2 + 55 + 35 2 + 5 3 3 2 + 65 + 35 2 + 5 3 

1 + 25 + 45 2 + 5 3 , ^ 1 + 2s + 4s 2 + s 3 

1 + 5s + 4s 2 + s 3 3 14 1 + 6s + 4s 2 + s 3 3 

1 + 2s + 4s 2 + 5 3 _ 1 + 2s + 4s 2 + J 3 

= 2 + 5s + 3s 2 + s 3 3 2 + 6s + 3s 2 + s 3 * 
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Figure 4. 


Theorem 11. IfG(s) is a stable transfer function, and<f> belongs to the sector [0, A:], 
then a sufficient condition for absolute stability is 

Re(- + g(joj)) > 0, for all u; G R 

rv 

(i.e j + g(s) is SPR). 

Combining this with our previous results we have the robust version of the Lur’e 
problem shown in Figure 5 below. 



Figure 5. 


Theorem 12. [21] Given the interval family G(s) of stable proper of stable proper 

systems and the family of sector bounded nonlinearities <j> belonging to the sector [0, k], 
a sufficient condition for absolute stability of the closed loop system is that k > 0 is 
any number such that 

k < oo, if inf inf Re(G(ju ;)) > 0 

Gjf (5) wEfi 


otherwise 

fc < 

infG*-(.) Re{G(juj)) 

where Gk(s) is the set of sixteen Kharitonov systems corresponding to G(^) . 
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This theorem may be generalized as follows. 

Theorem 13. Given the interval family G(«s) of proper systems stabilized by a fixed 
system F(s) 

inf inf Re[F(ju)( 1 + G(ju>)F(jiv)) 2 ] = 

G(s) 

inf inf Re[F{ju)){\ + G(ju)F(juj)) 1 ] 

In the last section we describe some frequency domain extremal properties of the 
extremal segments. 

8. EXTREMAL FREQUENCY DOMAIN PROPERTIES OF EXTREMAL 
SEGMENTS 

Consider again the feedback system shown in Figure 1. Since the extremal subset 
characterizes the robust stability of the interval system of Figure 1 it is natural to 
expect that these subsets also bound the Nyquist and Bode bands of interval systems. 
This is indeed the case and in this section we present recent results from [8] in this 
direction. This result was also independently reported in [22]. We expect these results 
to play a very significant role in synthesis and design issues. 

We shall give a quick summary of these results. 


8.1. Nyquist Envelopes 

Referring to the control system in Figure 1 we calculate the following transfer func- 
tions of interest in analysis and design problems: 


vW 

"( s ) 


<?« 


“M 

e(s) 


H>) 


H<) 


T y « 


^ = G( s) F W 

e(s) 

y( s ) _ GjsjFjs) 

• r(s) 1 + G( S )F( 5 ) 

f(f) 1 

' r{s) 1 + G(s)F( 5 ) 

«(*) _ no 

• r(s) 1 + G(s)F(sY 


(34) 


(35) 


As (j(^) ranges over the interval uncertainty set G(s) (equivalently, ( N(s ), D(s)) 
ranges over Af (s)xU(s)) the transfer functions T°(s), T y (s), T u (s), T e (s) range over 
corresponding uncertainty sets T°(s), T y (s), T u (s), and T*(s), respectively. In other 
words 


T °(s) := (G( a )F(«) : G(s) 6 G (s)} 
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( 36 ) 


T>W 

G(s)F (.) 

' ' 1 + G(*)F(,) 

T"W 

" X l + G(s)F(s) 

*■(<) 

■= { 1 

X l + G(s)F(a) 


GW € G(„)} 
GW € GW> 
G(t) € GW}- 


The extremal subsets of the transfer function sets (36) are also introduced: 


T£(-) 

TSW 


TSW 

T|(*) 


:= {G(s)F(s) : G(s) G G E (s)} 

r g(g)£(f) 

• l l + G(j)F(j) 


{ 


F(s) 


: G(s) G G E (s)} 
: G(s) G Ge(-s)} 


1 + G(s)F(s) 

b + GWFW : GW€GeW} 


(37) 


In frequency domain analysis and design problems the complex plane image of 
each of the above sets evaluated at s = jui plays an important role. We denote each of 
these two dimensional sets in the complex plane by replacing s in the corresponding 
argument by lj. Thus, for example, 


t eM '•= { t e( 5 ) : 5 = J"} (38) 


The Nyquist plot of a set of functions (or polynomials) T(s) is denoted by T: 


T := Uq^ooTH 


(39) 


The boundary of a set S is denoted dS. 

We shall give the main results here without proof. The proofs are given in Keel, 
Shaw and Bhattacharyya [23] and also independently by Tesi and Vicino [22]. 


Theorem 14. [8] [2] For every u> > 0, 

SG(u>) C G e (w) 

3T°(u>) C T|H 
dT v {w) C T£(w) 

0T u H C T£(w) 

3T» C T e E (u;) 

This result shows that at every ui > 0 the image set of each transfer function in 
(37) is bounded by the corresponding image set of the extremal segments. 

The next result deals with the Nyquist plots of each of the transfer functions in 

(37). 
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Theorem 15. [8] [2] The Nyquist plots of each of the transfer function sets T°(s), 

T y (s), T u (.s), and T e (s) are bounded by their corresponding extremal subsets: 

dT° C T E 
3T y C T v e 

dT u c T u 

dT e C T e E 

8.2. Bode Envelopes 

For any function say, T(s ) let //j(o>) := \T(ju>)\ and 0x(o;) := lT(ju>) denote the 
magnitude and phase evaluated at s = ju>. If T($) denotes a set of functions we 
let the extremal values of magnitude and phase at a given frequency be defined as 
follows: 


V- t(«*0 : = 

sup \T(juj)\ 

T(jw) 


p T (u>) := 

inf |T(w)|. 

T(jw) 

(40) 

Similarly 

^t(^) : = 

sup ZT(jaj) 
T(ju) 


; = 

inf IT(ju). 
T(ju) 

(41) 

Suppose that G(«s) is an interval family. To compute 


£ gM, 

£g( w ) 

(42) 

and 

4>g(v), 


(43) 

the following two lemmas are necessary. 


Lemma 5. Let A. be a closed polygon in the complex plane, and “a” be an arbitrary 
point in A. Let V \ be the set of vertices and E ^ be the set of edges of A ■ Then the 
following statements are true. 

1) max |a| = max |a| 

A V A 


2) min ja| = min |aj 

A E a 

Lemma 6, Let A and B be disjoint closed polygons in the complex plane, and “ a ” 
and <e b” be arbitrary points on A and B, respectively. Let Vji and Vs be the sets of 
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vertices and let E A and E B be the sets of edges of A and B, respectively. Then the 
following statements are true. 

1) max { la — Lb} = max {La — lb} 

} AxB 1 V A xV B 1 


min {la lb} = min {la — lb} 

AxB V A xV B 


Proofs of the above two lemmas are obvious from geometric considerations illus- 
trated in [2]. 

Let ff(w) denote the complex plane image of the set of polynomials N(s) 6 
evaluated at s = jw. Similar definitions hold for and Sj>(co). Ff(uj) is 

bounded by the set of Kharitonov segments Similarly, V{u) is bounded by 

the set Sp( u;). These facts along with Lemmas 5 and 6 lead to the following results. 
Before we state Theorem 16, let us define the following sets. 


GM 


Ge(w) 


G e (u;) 


(G(jw) 


N(juj) 


D(juj) 

N(jw) € JV(w)» D(joj) e V(u)} 

N(j") , 


{G{juj) 

N( 

{F(ju,) = 


D(ju) 

N(ju > ) € lC/s(u)), D(ju}) e St>(uj)} 

N(Ju) | 

D(ju) 

N(ju > ) € D(jio) G Kx>{v)}- 


(44) 


(45) 

(46) 


Theorem 16. For every frequency u> > 0, 

p G (uj) = £g e M 

£gM = £g b H 

Let us also define the set of systems constructed from Kharitonov vertices as 
follows: 

GkH := {F(juj) = | N{jv) 6 £«(<*), D(ju>) € K v {u>)}. (47) 

Theorem 17. For every frequency u> > 0, 

4>gH = <i>G K {u) 

4g( w ) = 

Using the above extremal properties it is possible to evaluate the Bode magnitude 
and phase bands of interval transfer functions. Let us consider the family of transfer 
functions 

T» = { 7» | F(»)G(.), G(«) € G(s) } . (48) 
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Since F(s ) is fixed, 


fi T o(uj) = \F{ju>)\ (ig{v) 

H T o(“) = \FU U )\ £g( u )' ( 49 ) 

Similarly, 

^T°( w ) “ iF(juj) + ^G^) 

<£ t » = iF(juj) + t G {u). (50) 

These relations are sufficient to construct the Bode magnitude and phase en- 
velopes. 

The Nyquist and Bode envelopes are important tools for solving analysis and 
design problems in robust parametric stability. In the next section, we show how the 
previous theory can be used to develop techniques to improve a given controller, by 
choosing an controller from a given set of stabilizing interval controllers, that provides 
optimal gain (or phase) margin to the closed loop system. 

9. DESIGN OF LINEAR INTERVAL CONTROL SYSTEMS 

In this section, we consider a nominal plant connected to an interval controller and 
give some design techniques for improving the closed loop gain and phase margins 
using the Nyquist envelope described in Theorem 15. From the results of the previous 
section it is clear that the main computational task is to determine the stability 
margin over the extremal segments. In the next section, we discuss the problem of 
determining optimal gain and phase margins over a single segment system. 


9.1. Segment System 


The typical extremal segment is of the form 


p(«) ; = { 


P o(a) 

Pi(s) + Ap 2 (s) 


AG [0,1]} 


(51) 


p(.) := |A€[°,1]} (52) 

Po{s) 

where are fixed polynomials. In this section, we develop simple techniques to 

compute the extremal gain and phase margins over a segment. We also determine 
the optimal value A*, equivalently p*(s), that produces the optimal gain (or phase) 
margin over the family p(«s). 

Let us consider the following segment system with 


<■ \\ Pofa) 

p JU> ’ Pi(jv) + Ap 2 (jw) 

The problem of computing the extremal gain and phase margins at the loop 
breaking point “m” over the single segment system is described as follows. Let us 
denote 

(Axfi) := {(A,w) | Lp{jw, A) = 180°, A 6 [0,1]} (53) 
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+ 


o 


m 


p(5,A) 


Figure 6. Segment System 


and 


Mp 

^ P 

<{> 

— P 

Then 

maximum gain margin over 

minimum gain margin over 

maximum phase margin over 
minimum phase margin over 


|p(j“>,*)l 

(54) 

|p(jw,A)| 

(55) 

Lp(ju, A) 

(56) 

Zp(ju;, A) . 

(57) 


pOO 

: P 

1 

(58) 

p(s) 

' P 

1 

Pp 

(59) 

POO 

: 6 

:= 4> p - 180° 

(60) 

p( s ) 

: 0 

:=4> -180°. 

—P 

(61) 


:= max 
(AxO) 

:= min 
(Axfi) 

: = max: 

(AxQ) 

:= min 

(AxO) 


Similar definitions can be made for the case of gain margins less than 1. 

As seen from eqs. (58) - (61), the problem of computing the extremal gain or 
phase margin over the segment system is two parameter optimization problem. This 
can be reduced to a simple one parameter problem as follows. Write 


Pi(ju) := Pifi( w ) + JPi/(^) 


Then 


p(j">*) 


Po(j^) 

pi(j^) + Ap 2 (jw) 


Poh(^) + jpoj(^) 

[PlilM + ^P2fl(<*>)] + i[Pl/( W ) + ^P2/(w)] 
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PorMpirM + Poi(<*QpiiQ «0 + A[pqr(o;)p2r(^) + Po/(^)p2f(^)] 

[pir(w) + Ap 2 fl(u;)] 2 + [pi/(^) + Ap 2 /(w)] 2 

' 

Re{p(ju/,A)} 

. Poj(^)pir(o>) ~ Por(^)pi/(^) + A[p 0/ (u;)p2fl(^) ~ Por(^)P2/M] 

J [pirM + Ap 2 fi(u;)] 2 + [pi/(o;) + Ap 2 /(u>)] 2 

v- 

Im{p(iw,A)} 

(62) 


In order to determine the gain margin, we set 

Zp(ju;, A) = 180° (63) 

which implies 

Im{p(ja;, A)} = 0. (64) 

Note that (64) will be satisfied when Zp(ju/, A) = 0° or 180°. We exclude frequencies 
u> for which Zp(jfa>, A) = 0°. From eqs.(64) and (62), we have 

Im{p(ju;, A)} = [po/(w)piji(u>) - Por(w)pu(c<;)] 

+A[po/(u;)p2R(c«;) - Por(<*0p2z(w)] 

= 0 


equivalently 


X(u) 


Por{u)Pii(v) ~ PojMpirM 


(65) 


(66) 


Po/(“>)P2r(“>) - Por(^)P2/(w)' 

From this representation, we can easily conclude that instead of searching both u; £ 
[0, oo) and A £ [0,1], searching only selected ranges of a ; that satisfy A € [0,1] is 
enough. Thus, we let 


A(tu) 


Por(^)Pi/(^) ~ Poj(^)Pir(^) 
Po/(^)P2i?(w) - Por(^;)p2/(^) 
— 0 or L 


(67) 


Without loss of generality, we have 

for A = 1 Por(v)Pii{v) ~ Poi(v)pir(u) - Poi(v)p 2 r(w) 
+Por(v)P2i(v) = o 

for A = 0 Pon(^)Pi/(^) - Poi(v)Pir(u) = 0 (68) 

The valid ranges of a; with respect to the condition A £ [0, 1] can be easily determined 
from the roots of the above two equations. Thus, the problem posed in eqs. (58) and 
(59) is reduced to selection of maximum and minimum magnitudes of A evaluated over 
the admissible ranges of u> determined from the roots of eq. (68). Furthermore, the 
optimal value A*, equivalently optimal values of parameters over the segment system, 
can also be easily determined by substituting uj* that corresponds to the maximum 
gain margin into eq. (66). 
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If the segment system is of the form in eq. (52), one can follow a similar procedure 
to determine the extremal margins and the corresponding optimal systems over the 
segment system. Similar procedures can also be applied for computing extremal phase 
margins over a single segment. This is easily derived by replacing the condition (63) 

by 

|p(ju>,A)| = 1 (69) 

9.2. Optimal Parameter Selection 

Applying the procedure described in the previous section to the entire set of segments 
systems, the extremal margins over the interval plant are determined. Consequently, 
the optimal system that produces the maximum gain or phase margin over an in- 
terval family is also determined. This procedure may be used to solve the following 
interesting problem. 

Suppose that a fixed system F(s) and a family of controllers G(s) are given, 
for which the closed loop system is stable. The objective is to select an optimal 
system G opt (.s) € IR n (s) so that the resulting closed loop system has the maximum 
possible gain margin or phase margin over the family G(s). Once such an optimal 
system is found the controller may be reset to the optimal parameter as the new 
nominal controller. At this point a new family of stabilizing interval controllers can 
be determined and the previous procedure of selecting the best controller repeated 
over the new box of parameters. The set of stabilizing interval controllers can be 
determined by many different methods; for example the locus introduced by Tsypkin 
and Polyak [24] may be used. This procedure described above can be repeated until 
1) improvement of the maximum margin in a given iteration is small or 2) the stability 
radius in the parameter space is small. Of course there is no guarantee that a globally 
optimum or even a satisfactory design will be achieved by this method. 

In the next section, an illustrative example is given. 


9.3. Illustrative Example 

Suppose 

F(s) := 

GW := 

where 


Ngjs) _ s 2 + 2s + 1 
D g (s ) s 4 + 2s z + 2 s 2 + s 

n x s + n 0 
d 2 s 2 + diS T do 


n 0 G [0.9, 1.1], mG [0.1, 0.2] 

d 0 G [1.9, 2.1], d x £ [1.8, 2.0], d 2 G [0.9, 1.0] 

We first check the stability of the family of closed loop systems with F(s) and 
G(s). This can be done by checking the stability of the corresponding extremal 
segment. 
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We have 

K^(s) = 0.1s + 0.9 K 2 n {s) = 0.2^ + 0.9 

kI(s) = 0.15 + 1.1 Kl(s) = 0.25 + 0.9 


Kfa) = s 2 + 1.8s + 1.9 Kl{a) = s 2 + 2s + 1.9 

= 0.9s 2 + 1.8a + 2.1 if3(s) = 0.9s 2 + 2s + 2.1 


and the corresponding segments Ge(s) are 

A*iM + (l-A)K>(«) Kj(.) 

K}(») AJfj(») + (1 - A )K’(s) 


where k = 1, 2, 3, 4 

(hj) £ {(1, 2), (1, 3), (2,4), (3,4)} 

and A € [0,1] Using the Segment Lemma [13], we verified that all the above extremal 
segments stabilize the closed loop. 

The Bode and Nyquist envelopes associated with T°(s), the forward transfer 
functions i^(s)G(s), are constructed by evaluating the following rational functions 
over o> € [0, oo) and A 6 [0, 1]: 

D(»)[AJfj(») + (1 - A)Xj(*)] 

and 

N(.){\Ki(s) + (1 - A)X'(.)] 

D(s)K$(s) 

This is shown in Figures 7,8 and 9. 

From these figures, the minimum and maximum gain margins are found to be: 

p = 1.1240. p = 1.7582. 

If we want to improve the gain margin of the system by selecting parameters and 
di beyond its previously given intervals, we can repeat the procedure as follows. With 
the controller designed in the previous part as the new nominal controller we can again 
construct an interval family of stabilizing controllers centered at the parameters of 
Gopt('S). This can be done by several methods. We used the stability locus introduced 
by Tsypkin and Polyak [24] which is shown in Figure 10. Note that in this case the 
locus shows the t 2 stability margin. 
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L.5 -1 

Figure 9. Nyquist envelope plot 



Figure 10. Stability locus (Tsypkin - Polyak) 

From this method, we have the parametric stability margin, 7 = 0.2128. Thus 
construct the stabilizing intervals around the nominal values and d{ as follows: 

7 7 

Ui - - < Tli < Tli + - 


di — ^ < d{ < di + -J- 




Consequently, we obtain 


_ + itj , 11 

< 4'W4'W<£ )> 

where 

n£ 1} G [0.7936, 1.0064], € [0.0839,0.2967] 

4° € [1.9936,2.2064], G [1.8936,2.1064], 4° G [0.7936,1.0064], 

Now the previous optimization procedure can be applied to this new family of interval 
controllers. This yields the result: 

AtfJM + (1 - X)KI( S ) 

G^s) = Ia=. 

0.2967s + 0.7936 
0. 7936<s 2 + 2.1064s + 2.2064 

and the maximum gain margin obtained from this first iteration is 

pW = 2.4295. 

The corresponding phase margin is 71.4506° (i.e., clockwise rotation of 71.4506 de- 
grees). 

By repeating the same procedure until the relative improvement of the gain mar- 
gin becomes small enough, we can obtain the “optimal” selection of the parameters 
for this problem. 

.3359s + .0398 

G o P t( s ) - 0398s 2 + 2.8602s + 2.9602* 

The maximum gain margin obtained from this second iteration is 

p* = 562.3651 

and the corresponding phase margin is 95.77° (i.e., clockwise rotation of 95.77 de- 
grees). Figure 11 shows the Nyquist plot of the corresponding optimal system for 
each iteration. 

10. MULTILINEAR INTERVAL SYSTEMS 

In following sections of the report, we consider systems in which the uncertain param- 
eters enter the characteristic polynomial affine multilinearly. As an example consider 
a feedback control system with a fixed compensator connected to a cascade of two 
interval plants as in the block diagram below: 

We have the following expression for the characteristic polynomial 

£(*) = F 1 (s)P 11 (s)P 12 (s) + F 2 (s)P 21 (s)P 22 (s) (70) 

with Pij(s) being interval polynomials i = 1,2 ,j = 1,2. The Theorem derived in [10] 
cannot deal with the robust stability of the family (70) because it contains products 
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Figure 1 1 . Nyquist plot for various controller parameteers. 



Figure 12. Unity Feedback System. 


of interval polynomials. Neither can the stability of the family (70) be checked by 
using the Edge Theorem of [12] since it is not a polytope. 

The above considerations motivate the problem of determining the Hurwitz sta- 
bility of the family of polynomials 

6(s) = F 1 (s)P n {s)P u (s)---P lri ( S ) + 

+Fm(s)P m i(s)P m2 (s) ■ ■ • P mr Js) (71) 

where Fi(s ) are fixed, and the polynomials Pij(s) are interval polynomials. The un- 
certainty set is therefore a box II in the space of these coefficients. The family (70) is 
of this type represents the special case of (71) when products of uncertain polynomi- 
als do not occur. Characteristic polynomials of the form (71) always occur in control 
systems containing several interconnected subsystems with uncertain parameters. We 
remark here that the vector of uncertain parameters, namely the set of coefficients 
of the polynomials Pij(s), enters into the characteristic polynomial coefficients affine 
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multilinearly. Since we assume that these parameters vary within prescribed inter- 
vals we refer to the family (71) as a multilinearly parametrized interval family. This 
form of the characteristic polynomial occurs in state space descriptions with interval 
matrices and also in matrix fraction description of multivariable systems when the 
matrix factors contain interval polynomials. 

We first introduce some notation. For any positive integer n let n denote the set 
of integers {1, 2, • • • n}. We consider the family of polynomials 

8{s) = F l (s)P u (s)P 12 (s)---P lri (s) + 

+ F m (s)P ml (s)P m2 (s) • ' ’ Pmrjs) (72) 


where Fi(s ) are fixed, and the polynomials Pij(s) are interval with i E m, j € r { . Let 
d°(Pij ) denote the degree of Pij(s) and let p l {j denote the coefficient of s l in Pij(s): 


PM = 4 + Ph* + • • • + Pi/ + - + 

Write p { . := (p°-, p}j, • • •pf/' 3 ' 1 ) for each i E m,j E r { and let 

P : =[£n>Pi2-- 


•P 

-mr ra J 


(73) 


(74) 


denote the vector of coefficients of the polynomials PM, i E rn,j E r i . Each such 
coefficient belongs to a given interval: 


p\j ^ [a-,-, 0 ij\ l = 0 , • - ■ , d°(P;;)> * G j G £;• ( 75 ) 

In the space of these coefficients we have the box II of uncertain parameters: 


n := {pIp-j € [a-j, 0^,1 = 0 > ”- > <f , (Py)»* £ £;}• (76) 

Each point p E II corresponds to a particular choice of the ordered set of polynomials 
Pij(s),i E m,j E r,-. We write 


£( 5 ,p) = f 1 {s)p u (s)pM---PiM + 

+ F m (s)P ml (s)P m2 (s) ■ • ■ P mrm (s) (77) 


to display the explicit dependence of £(s) on p. For a given fixed set of polynomials 
[Fi(s), F 2 {s) * • • F m (s)] := F_ let A denote the family of polynomials generated by 
the map F_ : II => 6 as in (72) and obtained by letting the parameter vector p, 
(equivalently, the polynomials P{j(s)) : range over the box II described in (76). In 
other words 

A := {S(s,p)\p e n}. (78) 

The four Kharitonov polynomials associated with the family of interval polynomials 
corresponding to Pij{s ) are 

KM = a° j + a} jS + (3y+fi jS * + --- 

KM = + <4-* + ■•• 

KM = + + + + 

KM = %+&*+<*?/ + <**/ + •■■> 
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and these definitions hold for each i E m,j £ Hi- Corresponding to each Pij(s) we 
define the four polynomial segments 

SJ := [*£(«), 1. s h -■= *SW1 

sfi ■■= (*£(.), JfJW), sj := jtJW), 

which we call the Kharitonov segments. These segments were introduced originally in 
[10]. A typical element of the segment Sjj, for example, is a polynomial, denoted by 
S}j(X,s) which is a convex combination of the form called the Kharitonov segments 

(1 - A )tf?.(s) + A Kh(s) := Sj( A, a), A € [0, 1], 

We now need to introduce some special subsets of II called the extremal manifolds. 
Fixing i — / we let 11/ C II denote the union of all the rj dimensional linear manifolds 
obtained by letting Pij(s),i / l range over the corresponding Kharitonov polynomials 
Kfj(s),k 6 4 ,j G Zi,i G m, and Pij(s), range over the Kharitonov segments Sjj, k G 
4 ,j G L t . The extremal manifolds Hi, l G m map into the corresponding multilinear 
surfaces A; C A, / G m in the space of coefficients of the polynomials £(s), under the 
previously defined mapping. More concretely, 


n, := 

- If If *( 2 ' 1 ) - . . If'( 2 > T *)( 


and 


Equivalently 
Finally, let 


{[^i J («) * • - K l \ • ■ • *?*(*)> ■ 

sg ,,,) P..*)sg , ' ,) p., *) • • • sg t " ,) (x 1 , •), , 

i(k, n) G 4, k G m, n G rk, A j G [0, 1], j G n } 

: = {<(-.) = .^"’w 

+ • ■ • + F,( S ).S; 1 (, ' 1 | (A 1 ,i).s/'‘' 2 | (A„ *).... .S^’t A,„ «) + 
| i(fc,n) G 4, fc G m, « G rjt, Aj G [0, l],j G n } 

A, := {£(a,p)|p G n,}. 


n E := [Jn, 


i=i 


denote the set of all linear extremal manifolds and let the corresponding set of mul- 
tilinear manifolds 

m 

A e := U A * = {^( a >P)IP € n E }- 

i=i 
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Because there is a one to one correspondence between the elements of IIe and of Ae 
we refer to both sets as extremal manifolds. Define the Kharitonov vertices , K(II), of 
II to be the subset of all vertices of II corresponding to the Kharitonov polynomials 
of the Pij(s). 

To illustrate the definition of the manifolds consider the special case of the family 
(70). In this case each manifold in II! is the union of polynomial vectors of the form: 

$,(*!, «), s&(Aj, .), KM, KMh (™) 

where X t G [0, 1], t G 2 and (z, j, k , l) range over (4 X 4 X 4 X 4). Similarly n 2 consists 
of the union of the polynomial vectors 

[^ll( 5 )j ^12( 3 )) *^2l(^l> s )i *^ 22 (^ 2 ? 5 )] (80) 

where (i>j, k,l) range over (4 x 4 x 4 x 4) and X t G [0,1], t G 2. The polynomial 
manifold consists of polynomials of the form: 

€ [0,l],t e 2, (81) 

and Ai consists of the union of such manifolds obtained by letting (z, j, l) range 
over (4 X 4 x 4 x 4). Similarly, the manifolds contained in A 2 are of the form: 

F 1 (j)/fS 1 («)Xi( < ) + F 2 (j)5j 1 (A 1> «)55 2 (A 2l 5) > A l <E [0,1], < <E 2, (82) 

where (z, j, A, /) range over (4 x 4 x 4 x 4). 

It is not difficult to see that the number of distinct manifolds in II/ in the most 
general case when all the Kharitonov polynomials associated with each polynomial 
Pij(s) are distinct, is 4( r i+ r 2+-+ r ”0, Since this holds true for each / G m the total 
number of extremal manifolds in IIe or Ae is 77x4^, R = + r 2 + • • • r m . With these 

preliminaries we are ready to state the main result of the next section. 

11. STABILITY OF MULTILINEAR MANIFOLDS 
11.1. The Multilinear Theorem 

In this section we give necessary and sufficient conditions for the Hurwitz stability of 
the family (71). Using the notation introduced in the last section we shall say that 
F_ stabilizes the family II if and only if each polynomial of the family A is Hurwitz 
stable. Similarly we shall say that F_ stabilizes IIe if an d only if every polynomial in 
Ae is Hurwitz stable. 

Theorem 18. (Multilinear Theorem) F_ stabilizes II if and only if F_ stabilizes 
IIe- 

The proof of this theorem is based on induction and may be found in [19]. 

Remark 1 . The assumption of independence of the perturbations can be easily re- 
laxed. The reader is is referred to [19] for the details. 
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11.2. Simple Determination of Stability of Two Dimensional Multilinear 
Manifolds 

In this section we consider the problem of checking the stability of a extremal manifold 
of dimension 2. This case will arise when = 2 in (71) and is interesting because it 
can be solved analytically. Consider therefore the following two dimensional manifold 

<5^(Ai , A 2 , -s) = Po( 5 )^1^2 + Pl( s )^l + P2( s )^2 + Pz{ s ), 

Ai € [0,1],*= 1,2. (83) 

Assuming the four vertices are stable, the manifold £(Aj, A 2 ,s) is unstable if and only 
if it has a jw root for a set of real values (A 1} A 2 ) G [0, l]x[0, 1]. To test for this we 
separate (83) into real and imaginary parts after substituting .s = ju, and set them 
equal to zero. This gives 

P0r( w )^1^2 + Plr( a ’)^l + P2r( UJ )^2 + VZri}* 3 ) = 0 

P0,(u>)A 1 A 2 + Plt(k>)Ai + P2i(w)A2 + P3i{u) = 0 (85) 

where 

Pfc(s)|.=io; :=pfcr(w) + jPki{w), for fc = 0, 1,2,3. 

From (84), we have 

[po r (u>)A 2 + pi r (u;)]A 1 + p 2r (^)A 2 + P3r(^) = 0 (86) 


and 


P2r(^)A2 + P3r(^) 

POr(w)A 2 + Plr(^) 


(87) 


Similarly, from (85), 

» P2»(<*QA 2 + P3»(^) / g g \ 

P0i( w )^2 + Pli^) 

Since A a = oo ^ [0, 1] we can without loss of generality, deal only with the case in 
which the denominators of (87) and (88) are nonzero. By equating (87) and (88) we 
have, 

[P2i(w)pOr(w) - Pot(w)p 2 r(w)]A 2 

+ [p 3i (uj)p 0r (u}) + p 2i (a>)p lr (u>) - p 0i (u!)p 3r (u;) - p li (u;)p 2 r (u;)]A 2 

+ [P3i(w)Plr(w) - Pli(w)P3r(^)] = 0 (89) 

From (89) and (87) we can solve for Ai(u;) and A 2 (u;) and verify if (A^a;), A 2 (u;)) 
intersects the set [0, l]x[0, 1] for some u;. If the intersection is empty the manifold is 
Hurwitz stable otherwise it is unstable. 
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12. EXTREMAL PARAMETRIC STABILITY PROPERTY: MULTI- 
LINEAR CASE 

We now consider the family A and the manifolds lie and Ae which occur in the 
Multilinear Theorem of the last section. As before let 

P ;= fell* £l2 PmrJ 

denote the n dimensional parameter vector consisting of the ordered set of coefficients 
of the polynomials Pij(s) and let p £ R n vary in the prescribed box II specified by 
the given upper and lower bounds: 

v[j e I = 0 , • • • , d 0 (Pij), i<Em,j e r { 

Let || • || denote any norm in R n and let V u denote the set of points u in R n for which 
6(S,u) is unstable or loses degree (relative to its generic degree over II. Let 

p ( p ) = Up - u ll p 

denote the radius of the stability ball (measured in the norm || • ||) and centered at 
the point p. This number serves as the stability margin associated with the point p. 

If the box II is stable we can associate a stability margin with each point in II. A 
natural question to ask then is: Is there a point in II which is closest to instability 
in the norm || • || and where is it? The answer to that question is provided in the 
following theorem. 

As before we define a mapping from II to the set of all positive real numbers: 

n n + \{o} 

p — ♦ p{ p) 

Our question stated in terms of functions is: Has the function p(p) a minimum and 
is there a precise point in II where it is reached? 

Theorem 19. (Extremal property of the stability manifolds) The function 

n ft + \{o} 

p — ♦ p ( p ) 

has a minimum which is reached at a point on the extremal manifolds IIe- 
The proof of this theorem may be found in [19] and omitted here. 

13. PARAMETRIC AND UNSTRUCTURED PERTURBATIONS: MUL- 
TILINEAR CASE 

In this section we will analyze the problem of robust stability in the presence of both 
parameter variations and unstructured perturbations modelled in the usual way as 
norm bounded perturbations. The subject of robust stability under mixed types of 
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perturbations is of current interest (see for example [25], [26], [20], [21] and [27]). We 
model this situation by considering a multilinear interval plant, namely one whose 
transfer function is a ratio of polynomials of the type that was introduced in 71. To 
be specific we will consider single-input, single-output, proper, stable systems with 
transfer function of the form 

s) = 

Here 



7 (*) = H 1 (s)L u (s)L 12 (s)-^L lri (s) + 

T (s)T m 2(s) * • “ 

where the polynomials Hi(s) are fixed and the polynomials L{j(s) are interval polyno- 
mials, that is their coefficients vary in a prescribed box A; the corresponding family 
of polynomials 7 (s) is denoted by I\ We suppose as before that 

6(s) = F 1 (s)P ll {s)P 12 (s)---P lri (s) + 

where the polynomials F t (s) are fixed, the polynomials Pij(s) are interval polyno- 
mials, with coefficients that vary in the prescribed box II and the resulting family 
of polynomials £(s) is denoted A. As in Section 2 we let p denote the vector of 
coefficients of the polynomials {P^(s)} and we similarly let 1 denote the vector of 
coefficients We also denote explicitly, the dependence of £(5) on p and of 

i/(s) on 1 by writing 5(s, p) and i/(«s,l) whenever necessary. It is assumed that the 
parameters p and 1 perturb independently. From these polynomial families we form 
the parametrized family of transfer functions 

Q = {ttA Ip € n, and 1 € A}. (90) 

<b(s,p) 

To display the dependence of a typical element g(s) of Q on 1 and p we write <7(5, p, 1). 
Introduce the Kharitonov polynomials and segments associated respectively with the 
Pij(s) and Lij(s) respectively. As in Section 2 these are used to generate the extremal 
subsets IIe of El and A E of A respectively. The Kharitonov extreme points of II and 
A are denoted respectively by K(II) and K(A). Finally, we denote the polynomial 
manifolds resulting from K(II), K(A), Ae and IIe as follows: 

r E = {7(-,l)|l € A E },r* - {7(M)|1 € K(A)} 

A e = {$(*,p)|p € n E }, Ajf = {A(js, p)|p € K(II)}. 

The main results in this section will deal with the calculation of the H stability 
margin for systems containing parameter uncertainty as defined above. In the follow- 
ing we will use the standard notation: C+ := {s 6 C : Re(s) > 0}, and H^C^.) 
will represent the space of functions f(s) that are bounded and analytic in C+ with 
the standard norm, 

WfWoo =sup|/(ju>)|. 

u£R 
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To determine the unstructured stability margin of the family Q we need to deter- 
mine the supremum of the norm of certain transfer functions over Q. Specifically 
we formulate the following problems: Let W{s) be a scalar stable weight, with a 
stable inverse, and write W(s) = . 

A) Consider the feedback configuration shown in Figure 13, Q is a stable family, and 
A P is any H ^ perturbation that satisfies ||AP|| < a. 




AP 


Figure 13. Multiplicative Perturbations 


B) Consider the feedback configuration shown in Figure 14, AP is any H ^ perturba- 
tion that satisfies || AP || < a, and C is a controller that simultaneously stabilizes 
every element in the set Q. 



Figure 14. Additive Perturbations 


The above problems are generalized versions of standard robust stability 
problems (see [28]) where a fixed plant is considered. The solution is accomplished 
once again by showing that the Poo norms in question attain their supremum value 
over a certain extremal set of transfer functions Q& C Q . This set is defined as follows: 

g E := {^4 1(1 £ k ( a )>p G n E) or C 1 e A E)P € K(n)}. 

We can now state the main result of this Section. 
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Theorem 20. (Extremal properties) 

A) sup ||Wflf||oo = sup ||Wj?||oo, 

g£Q 9 £Ge 

B) sup \\WC(1 + gC)- 1 Hoc = sup \\WC(l + gC)-^. 

g£G QtG e 

Corollary 2. (Unstructured Margins) 

1) The configuration of Figure IS will be stable if and only if a satisfies 

1 . 

— II II * a o- 

su P fl ef? B M°° 

2) The configuration of Figure 14 will be stable if and only if a satisfies 

1 = , 

a “™p oefe ||C(i + 9 c)-'|U' ■=' 

The proof of this theorem is similar to that used in [20] and details are omitted here. 
The idea behind this approach is to replace the question of finding an upper bound of 
the Hqq norm of a transfer function by an equivalent question concerning the stability 
of a certain parametrized family of polynomials, for which the results of the previous 
sections apply. For this purpose we need the following lemma [20] which gives a 
characterization of proper rational functions g(s) which are in /f oc ( C+) and which 
satisfy WgW^ < 1. 

Lemma 7. Let h(s ) = n(s)/d(s) be a proper (real or complex) rational function in 
# 00 ( 0 +), with deg(d(s)) = q, then ||&||oo < 1 if and only if 

al) | n q \ < \d q \, 

bl)d(s) + e^ e n(s) is Hurwitz for all 9 in [0,27t). 

Remark 2. The quantities a* and a* serve as unstructured stability margins for 
the respective open and closed loop parametrized systems treated in Problems I and 
II. 

14. PARAMETRIC AND NONLINEAR PERTURBATIONS: MULTI- 
LINEAR CASE 

Another stability robustness problem that involves structured and unstructured per- 
turbations is the classical Lurie problem of nonlinear control theory. This problem 
considers a fixed linear time invariant system subjected to perturbations in the form 
of nonlinear feedback gains contained in a prescribed sector. In [21] a robust version 
of the Lurie problem was treated. In this formulation the fixed linear system is re- 
placed by the more realistic model of a parametrized family of plants. The “nonlinear 
stability margin” of the system can be determined by finding the infimum, over the 
parametrized family, of such stability sectors. From standard results on the Lurie 
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problem, the size of such a sector can be determined by finding finding the infimum 
of the real part of g(juj) as g ranges over the parametrized family. In [21] it was 
shown how the strict positive real (SPR) property for a stable family of interval sys- 
tems can be determined from a set of sixteen plants called the Kharitonov systems. 
In addition, in the presence of a fixed controller that stabilizes an entire family of 
interval systems, the SPR property for the family of transfer functions (7(1 + gC)~ l 
is determined from a set of 32 one parameter family of systems. Here we consider 
the more general situation where the parametrized family considered is the family 
Q defined in the previous section. Using the extremal properties established in the 
last section and the proof developed in [21], it is possible to establish the following 
theorem. The proof is omitted as it is very similar to that of the last section. 

Theorem 21. (Extremal properties) 

1 ) Let Q be the multilinear family defined above, and assume that Q is stable then 

inf inf Re(W(ju>)g{ju))) = inf inf Re(W(jw)g(ju>)). 

gtzG u£R g€G e <*>£R 

2) If C is a controller that stabilizes the entire family Q , then 

inf inf Re(WC( 1 + gCy^jcv)) = 

g£G u>£R 

inf inf Re(WC( 1 + gC)~ l {ju)). 

g£G e w£.R 

15. CONCLUDING REMARKS 

The summary of results presented here form the beginnings of a complete theory 
of interval control systems. We expect such a theory to develop over the next few 
years. We expect such a theory to impact on the design of control systems filters and 
communication systems. 
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